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ABSTRACT 

The anti-correlations between the equivalent widths of emission lines and the continuum luminosity 
in AGNs, known as the Baldwin effect are well established for broad lines, but are less well studied for 
narrow lines. In this paper we explore the Baldwin effect of narrow emission lines over a wide range of 
ionization levels and critical densities using a large sample of broad-line, radio-quiet AGNs taken from 
Sloan Digital Sky Survey (SDSS) Data Release 4. These typel AGNs span three orders of magnitude 
in continuum luminosity. Wc show that most narrow lines show a similar Baldwin effect slope of 
about -0.2 while the significant deviations of the slopes for [Nil] A6583, [Oil] A3727, [Ncv] A3425 
, and the narrow component of Ha can be explained by the influence of mctallicity, star-formation 
contamination and possibly by difference in the shape of the UV-optical continuum. The slopes do 
not show any correlation with either the ionization potential or the critical density. We show that a 
combination of 50% variations in continuum near 5100A and a log-normal distribution of observed 
luminosity can naturally reproduce a constant Baldwin effect slope of -0.2 for all narrow lines. The 
variations of the continuum could be due to variability, intrinsic anisotropic emission, or an inclination 
effect. 

Subject headings: galaxies: active-galaxics:Seyfcrt-(galaxics:) quasars: emission lines 



1. INTRODUCTION 

The anti-correlation between the equivalent widths of 
broad emission lines and AGN luminosity (the "Baldwin 
effect" ; hereinafter BE) was first discovered by Baldwin 
(1977) for the ClV A1549 broad emission line in high 
redshift AGNs. It was initially hoped to be able to use 
the effect to calibrate the AGN luminosity to be able 
to use AGNs as cosmological standard candles (Bald- 
win et al. 1978), but the large dispersion of this rela- 
tionship rendered this impossible (Baldwin ct al. 1989; 
Zamorani et al. 1992). The BE is now well established 
for nearly all broad emission lines and the slope of the 
BE steepens with increasing ionization potential (Zheng 
& Malkan 1993; Dietrich et al. 2002). Several mech- 
anisms have been proposed to explain this effect (See 
Shields 2007 for a review). Among these perhaps the 
most widely accepted one is that the ionizing continuum 
softens with increasing luminosity so there are relatively 
fewer ionizing photons for broad emission line formation 
in high-luminosity AGNs. This model can reproduce the 
ionization energy-BE slope relationship fairly well (Ko- 
rista et al. 1998) and the assumption has observational 
support (Binette et al. 1989; Zheng & Malkan 1993; 
Wang & Lu 1998; Korista et al. 1998). Some theoretical 
models (Netzer 1985, 1987; Netzer, Laor, & Gondhalekar 
1992; Wandel 1999a, b) could produce a softer ionizing 
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spectrum in high-luminosity sources, but the standard 
thin disc model (Shakura & Sunyaev 1973) they adopt 
suffers from many problems (Antonucci 2002; Gaskell & 
Klimek 2003; Gaskell 2008; Lawrence 2012; Antonucci 
2012). In this sense, the explanation of broad line BE is 
still elusive and controversial. An outlier of the ioniza- 
tion energy-BE slope relationship is Nv A1240, which 
has an ionization energy of 97.7eV but shows no BE at 
all (Dietrich et al. 2002). It has been proposed that 
this can be explained by a dependence of metallicity on 
AGN luminosity, which in turn is a combination of Ed- 
dington ratio (L/LEdd) and black hole mass (Korista et 
al. 1998; Hamann & Ferland 1993, 1999; Dietrich et 
al. 1999; Dietrich & Wilhclm-Erkens 2000). In local 
galaxies, gas metallicity correlates well with the mass of 
galaxies (Tremonti ct al. 2004), and more massive galax- 
ies have larger black hole masses (-Mbh) according to 
the Mbh-M* relationship. For Typel AGN population, 
the brighter AGNs would have higher A/bh on average 
(Kollmeier et al. 2006; Steinhardt & Elvis 2010a, b; Lusso 
et al. 2012) and thus higher metallicities. Recent studies 
of the BE for C IV, Mg II and Fe II, however, challenge this 
picture by showing that the BE might instead be driven 
by the EWs correlation with Eddington ratio (Baskin 
& Laor 2004; Bachcv ct al. 2004; Warner, Hamann & 
Dietrich 2004; Zhou ct al. 2006; Dong ct al. 2009a,b) or 
with Mbh (Netzer, Laor & Gondhalekar 1992; Wandel, 
Peterson & Malkan 1999; Shields 2007; Kovacevic et al. 
2011). 

The BE of narrow emission lines is much less well 
studied, and some results are still controversial. Stcincr 
(1981) discovered a strong BE for [Olll] A5007 (espe- 
cially in AGNS with strong optical Fe II) and Wills et al. 
(1993) found the narrow lines in high-luminosity AGNs 
are very weak. Meanwhile, Croom et al. (2002) found 
a significant BE for [Ne V] A3425 and [O II] A3727 but 
obtained a null result for [Nc III] A3870 and [O III] A5007 
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using the 2dF sample. Subsequent work shows, however, 
that [O III] A5007 docs show a BE (Dietrich et al. 2002, 
Netzer ct al. 2004). Honig ct al. (2008) and Keremed- 
jiev et al. (2009), used Spitzer data to find that the BEs 
of different mid-IR narrow emission lines have a nearly 
constant slope. Several mechanisms have been proposed 
to explain the NLR BE. NLRs follow a size-luminosity 
relationship (Schmitt et al. 2003b; Bennert et al. 2002; 
Greene et al. 2011). The size of the NLR may grow 
beyond the size of the host galaxy in high-luminosity 
AGNs, thus turning from ionization-bounded to matter- 
bounded, and so producing the BE (Croom et al. 2002). 
As the EW is proportional to the covering factor (CF) 
and it is found that the CF contributes to much of the 
variance of the EW (Baskin & Laor 2005), a luminosity- 
dependent CF is also a possible cause of the BE (Shields 
et al. 1995; Stern & Laor 2012b; Stern & Laor 2012c). 
It has recently been proposed that the EW of narrow 
emission lines is dependent on the inclination to the ac- 
cretion disc (Risaliti et al. 2011), at least in the highest 
EW sources. In principle, this is a potential cause of the 
BE too. 

To make progress in understanding the NLR BE, we 
need to measure narrow emission lines of a wide range of 
ionization potentials and critical densities, and determine 
the nature of their BEs more accurately. In this paper, 
we use a well-defined sample drawn from SDSS DR4 and 
employ the technique of composite spectra to investigate 
the BE for prominent narrow emission lines in the optical 
band and to study the origin of the NLR BE. We use a 
cosmology with Ho = 70 kms -1 Mpc -1 , fi m = 0.3, and 
SIa = 0.7 throughout this paper. 

2. SAMPLE AND MEASUREMENTS 
2.1. Sample 

We need accurate emission line and continuum mea- 
surements to ensure reliable determinations of equivalent 
widths and continuum luminosities. From the spectral 
data set of the Sloan Digital Sky Survey Fourth Data Re- 
lease (Adelman-McCarthy et al. 2006), we have selected 
4178 Seyfert 1 galaxies and quasars (i.e., type 1 AGNs) 
as described in Dong et al. (2011). We apply a redshift 
cutoff of 0.8 so that the redshift of the spectrum can 
be accurately determined using [O III] A5007. To ensure 
high-quality spectra we require a median signal-to-noise 
ratio (S/N) of > 10 per pixel in the optical. To mini- 
mize the host-galaxy contamination (see the Appendix 
of Dong et al. 2011), we restrict the weak stellar ab- 
sorption features, such that the rest-frame EWs of CaK 
(3934 A), CaH + He (3970 A), and R5 (4102 A) absorp- 
tion features are undetected at < 2 a significance. This 
criterion ensures that the host contamination is less than 
10% around 4200A(Dong et al. 2011). One may note 
that Dong et al. (2011) do not consider very young stel- 
lar population (i.e., emission line galaxies). But previous 
analysis have shown that in massive galaxies, the optical 
continuum is not dominated by very young stellar pop- 
ulation although the UV continuum may be (Schawinski 
et al. 2007). After removing duplications and sources 
with too many bad pixels in the H/3 + [O III] region, we 
obtain 4178 type 1 AGNs. Including radio-loud AGNs 
may influence the measurement of EW of emission thus 
produce false effect for two reasons: firstly, the jet may 



interact with the ISM to enhance the narrow line emis- 
sion (Labiano et al. 2007) and secondly, it might also 
enhance the continuum through beaming if the jet points 
close to our line of sight. By matching with the FIRST 
catalog (Becker et al. 1995) using the method of Lu ct 
al. (2010), we reject 499 radio-loud AGNs so that our 
final sample consists of 3677 sources. 

2.2. Spectral fitting and measurements 

We give a brief description of our spectrum fitting pro- 
cess here; the details can be found in Dong ct al. (2011). 
To model the spectrum, we used a code based on the 
MPFIT package (Markwardt 2009) and fit the AGN fea- 
tureless continuum, the Fell multiplets, and other emis- 
sion lines simultaneously. The AGN continuum is repre- 
sented locally by a power-law, for the region of 4200-5600 
A and for the Ha region (if present). The Fell template 
by Veron-Cetty et al. (2004) that we use, is constructed 
using the identification and measurement of Fe II lines in 
I Zw 1. It has two separate sets of templates in analyt- 
ical forms, one for the broad-line system and the other 
for the narrow-line system. Within each system, the rel- 
ative velocity shifts and relative strength are assumed to 
be the same as those in I Zw 1. Broad Fell lines share 
the same profile as broad H/3, while each narrow Fell 
lines is modeled with a Gaussian. During the fitting, 
the normalization and redshift of each system are taken 
as free parameters. The broad Balmer lines are fitted 
with as many Gaussians as is statistically justified. All 
narrow emission lines, except for the [Olll] AA4959, 5007 
doublet lines, are fitted with a single Gaussian. Each 
line of the [Olll] doublet is modeled with two Gaus- 
sians, one accounting for the line core and the other for 
a possible blue wing, as seen in many objects. Since 
the sources in our sample do not suffer from significant 
host galaxy contamination, we do not apply starlight cor- 
rections to individual spectra. For each source, we use 
AL A (5100) and FWHM of broad H/3 line to obtain the 
A/bh using the virial mass estimates by the Dibai method 
(Dibai 1977) using the formalism of Wang et al. (2009). 
The typical statistical scatter about the A/bh obtained 
by reverberation-mapping is about 0.4 dex, and it may 
also be subject to more systematic errors. (Krolik 2001; 
Collin et al. 2006; Shen et al. 2008; Fine et al. 2008; 
Marconi et al. 2008; Denney et al. 2009; Rafiee & Hall 
2011a; Steinhardt 2011). The error of L/LEdd is of sim- 
ilar magnitude as that of A/bh- 

2.3. Composite Spectra Generating and Fitting 

A convenient way to explore the correlations between 
EWs and other parameters like luminosity is to make 
composite spectra for different parameter bins. We nor- 
malize individual spectrum to the mean flux around 
4200A and then construct the geometric composite as 
Vanden Berk et al. (2001). To get an accurate contin- 
uum measurement, we need to subtract the broad emis- 
sion lines, especially, Ha and H/3 and Fell emission from 
the spectrum. For the broad component of H/3 and Fc II 
we subtract them from the original spectrum; but for the 
broad component of Ha, we leave it un-subtractcd when 
making composite spectra. This is because the broad Ha 
is highly blended with [Nil] A6583 and [Nil] A6548, so 
deblending in individual spectrum is not reliable while 
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deconvolving the blends in the high SN composite spec- 
trum is easier. The fitting algorithm used to model the 
composite spectrum is the same as described in Section 
2.2. More specifically, we fit the blends using one gaus- 
sian for the narrow component of Ha, [Nil] A6583 and 
[Nil] A6548, two gaussians for [Sll] AA6717, 673f and 
three for Ha 6 . The broad-line-subtracted spectra are 
shown in Fig. 1. The embedded panels show the fitting 
result for the Ha fc + [Nil] blends. 

3. RESULTS 

3.1. The Baldwin effect for different lines 

First, we want to explore if the BE exists in the promi- 
nent narrow lines [NeV] A3425 , [Oil] A3727, [Nelll] 
A3870 , H/3 n , [Om] A5007, [Ol]A6300, Ha", [Nil] A6583 
and [Sll] AA6717, 6731. Measuring narrow lines in in- 
dividual spectrum are subject to the S/N limit and the 
EW-L5100 relation of all our lines shows a large dispersion 
(typical 0.2 dex, this can be seen clearly in the middle 
panel of Fig. 4), so we turn to composite spectra for reli- 
able measurements. We divide our sample into intervals 
of 0.3 dex in 5100A luminosity starting at log -L5100 = 
43.2 to 45.9 ergs~ 1 . and make a composite spectrum in 
each luminosity bin as described in Section 2.3. From 
Fig.l we can see two clear and remarkable results: 

• With increasing luminosity, the narrow lines van- 
ish. 

• The slope of the observed continuum become bluer 
with increasing luminosity. 

To see the dependencies of the EWs on the luminos- 
ity more clearly, in Fig. 2 we plot the log EW that de- 
rived from composite spectra for each luminosity bin for 
each of the lines listed above against log £5100 in or- 
der to get more qualitative results. For each line, we 
show a weighted linear regression of the logarithm of 
the EW on the logarithm of the luminosity and we list 
all the BE slopes of the narrow emission lines together 
with their ionization energies and critical densities in 
Table. 1. When fitting the relationship, we add an er- 
ror to each data point that would help to reduce the 
normalized \ 2 to ~1. These added errors are the same 
for individual narrow line, and they account for the po- 
tential other systematic errors. The BE slopes of the 
high-ionization lines [Nelll] A3870 and [NeV] A3425 are 
-0.26 ± 0.02, -0.31 ± 0.015 respectively which are con- 
sistent with the results of Keremedjiev et al. (2009) who 
obtained -0.22 ± 0.06, and -0.19 ± 0.06 respectively for 
the Nelll A 15.56^m, and NeV A 14.32/im lines. Our 
[Om] A5007 BE slope of -0.21 ± 0.016 is similar to the 
values of ~ —0.2 found by Kovacevic et al. (2011) and 
steeper than the —0.1 ± 0.02 found by Dietrich et al. 
(2002). The difference might be due to different sam- 
ple we employ. The error of slopes we give here include 
only the statistic error of the fitting, but do not include 
the measuring error and intrinsic dispersion of the EW, 
which may reach 0.2 dex typically. 

In addition to finding the BE in high-ionization lines, 
we also find the BE in the low-ionization lines of 
[Ol]A6300 [S II] AA6717, 6731 and [Oil] A3727. The nar- 
row recombination line: H/3™ also shows a BE with a 
similar slope as for the forbidden lines. For compari- 
son, the broad H/3 line shows an inverse BE with a slope 



of 0.16 for our sample. This is consistent with Greene 
& Ho (2005) and Croom et al. (2002) who reported 
EW{H/3 b ) (x L%$ and EW{H/3 b ) oc L%$ respectively. 
We can see in Fig. 2 that the EW of the narrow compo- 
nent of Ha increases with luminosity too. These indicate 
that the behaviors of broad lines and narrow lines are dif- 
ferent. 

The correlation coefficients between EW([OlIl]) and 
AL A (5100), z, M bh , L/L Ed d are -0.12, -0.05, -0.01, - 
0.21, respectively. These are similar to the results in 
Zhang et al. (2011). For our flux-limited sample, the 
rcdshift, luminosity, black hole mass, and Eddington ra- 
tio are correlated with each other. The correlation co- 
efficients between AL>(5100) and z, Mbh, L/Lsdd are 
0.81, 0.67, 0.15 respectively. After controlling for z or 
L/LEdd in a partial correlation analysis, EW([0 III]) cor- 
relates with AL A (5100) with r s =-0.14 and -0.16. This 
indicates that the BE of narrow lines is not a secondary 
effect of correlations between the EW and redshift or 
L/LEdd but an independent phenomenon. We note that 
as the correlation strength is not strong, one possibility 
is that there are other factors that regulate the EW as 
discussed in Zhang et al. (2011). Also, measurement 
error in EW([OlIl]) may act to smear the correlation. 
The measurement error could be introduced by limited 
S/N and fitting process. In spite of these uncertainties, 
the BE of narrow lines docs exist and may shed light on 
physical process in AGNs as explored in detail below. 

3.2. The ionization energy - BE slope relationship: no 
correlation 

A dependence of the BE slope on ionization energy 
is found for broad emission lines in AGNs (see intro- 
duction) and this is the most compelling evidence for a 
luminosity-dependent ionizing spectrum. The NLR lies 
far from the nucleus, has a complex geometry and may 
contain dust (Netzer & Laor 1993; Tomono et al. 2001; 
Radomski et al. 2003; Schweitzer et al. 2008). This 
makes the response of the narrow-line flux to changes 
in the SED complicated. We plot the slope of BE 
against the ionization energies of the narrow lines in the 
left panel of Fig. 3. We also plot the Keremedjiev et 
al. (2009) data for [Siv] 10.5lAtm:-0.29 ± 0.05, [Nell] 
12.81/mi:-0.25±0.06, [Nelll] 15.56/zm:-0.22 ± 0.06 and 
[NeV] 14.32^m:— 0.19 ± 0.06 in purple crosses and the 
broad-line BE from Dietrich et al. (2002) in blue rect- 
angles for comparison. We can see that the narrow-line 
BE slopes do not correlate with the ionization-cnergy ( 
Pnuii = 0.78 meaning the two-sided probability that a 
correlation is not present is 78%.) but cluster around 
-0.2 with a dispersion of ±0.1. [Nil] A6583, Ha" and 
[O II] A3727 have slopes of -0.10 ± 0.014, -0.29 ± 0.033, 
and — 0.37± 0.011 so they deviate from -0.2 significantly. 
Possible reasons for this are discussed in Section 4.1. 

3.3. Critical Density - Slope relationship: No 
correlation 

It is well known that the NLR is stratified that the 
high-ionization lines rise from the inner part of the NLR 
while low-ionization lines rise further away (see, for ex- 
ample, Veilleux et al. 1991; Robinson et al. 1994; Ben- 
nert et al. 2006a, b; Kraemer et al. 2009). The central 
electron temperature, density, and ionization parameter 
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are, in general, higher in Seyfert Is than in Seyfert 2s 
(Gaskell 1984; Schmitt 1998; Bennert et al. 2006b). The 
lines of lower critical density are not as strong near the 
nuclei as the higher critical density lines. A possible ex- 
planation of this is that the change in critical density 
marks a transition between difference types of clouds. 
To see if this is a factor in the NLR BE we therefore 
plot the EWs against critical density in the right panel 
of Fig. 2 using the same symbols of panel (a). We again 
fail to find any correlation in this plot (P nu ii = 0.82 using 
a Spearman rank correlation analysis). 

4. DISCUSSIONS 

Using our sample of 3677 radio-quiet AGNs from SDSS 
DR4 that have little contamination by host galaxy, we 
find that, in contrast with the BE of broad emission lines, 
the narrow lines show a nearly constant BE slope irre- 
spective of ionization potential or critical density. The 
[O II] A3727, [N II] A6583 and Ha™ lines show significantly 
different BE slopes from other lines. 

4.1. Deviation from constant BE slope 

Although the slope of different lines cluster around - 
0.2, we do find some lines that show large deviations. 

Firstly, [Nil] A6583 show a much flatter slope (-0.10± 
0.014) than other lines (P < 5% in two-sided KS test, 
meaning if the two samples are drawn from the same dis- 
tribution, we expect to see a difference as large or larger 
than we see here only 5 in 100 times.). This flattening 
was also found for N V A1240, which has an ionization po- 
tential of 77.7 eV but no BE. One plausible explanation 
of a weaker or absent BE for the nitrogen lines is that 
the increase of metallicity in brighter AGN will enhance 
nitrogen abundance relative to carbon, oxygen etc., be- 
cause nitrogen is a secondary element whose abundance 
scale as Z 2 . This trend could compensate for the de- 
crease of EW in higher-luminosity sources thus making 
the BE slope flatter. This explanation is supported by a 
number of recent studies. The metallicities of the BLR 
(Hamann & Ferland 1993; Nagao et al. 2006a; Juarez 
et al. 2009) and NLR (e.g., Nagao et al. 2006b; Mat- 
suoka et al. 2009) are both found to correlate with the 
luminosity of the AGN. So in principal, the Z/vl.r-L rela- 
tionship could produce the flattening of the BE of [N II] 
A6583 we observe here. However, this explanation has 
a major problem because Niv] and Nlll], which would 
be expected to deviate from the ionization energy slope 
relationship for the same reason, lie on it (Dietrich et 
al. 2002). Despite this serious problem to be resolved, 
Occam's razor suggests that an enhanced abundance in 
high-luminosity AGNs is the simplest and most plausible 
explanation of the flatter BE slope of [Nil] A6583. 

The second significant deviation from a BE slope of 
-0.2 is the steep slopes of Ha™ and [Oil] A3727. Both of 
these are star formation (SF) indicators in star-forming 
galaxies (Kenicutt et al. 1998, Ho 2005) and their EW 
may reach several hundred A in starburst galaxies. So 
even though the AGNs we selected show no absorption 
feature in the continuum, the emission lines are still pos- 
sible subject to SF contamination because of their large 
EW. In low-luminosity sources whose continuum is low, 
the emission line from SF may contribute significantly to 
the total line flux, and thus enhance the EW. In type-2 
AGNs, the SF contribution of Ha" is estimated to be 



more than 60% (Brinchmann ct al. 2004). The Ha 
emitting region of typc-1 AGNs, however, may have a 
higher fraction of the emission originating with the AGN 
(Zhang et al. 2008) but still be heavily influenced by SF. 
Furthermore, a variety of studies have established a cor- 
relation between the strength of AGN activity and star 
formation in the local universe ( e.g., Rowan- Robinson 
1995; Croom et al. 2002; Netzer et al. 2007, 2009; Shao 
et al. 2010). The steepest correlation is Lsf cx L°j^ n 

(Netzer 2009) meaning Lsf/Lagn °c -^agjv> a BE slope 
steeper than -0.2. So it makes sense that the SF in host 
galaxies would steepen the BE slopes of Ha™ and [O II] 
A3727. However, the [O II] A3727 and Ha™ BEs cannot 
be attributed entirely to SF-AGN relationship. For a 
plausible range of ionization parameters, densities, and 
ionizing spectra, the intensity of [O II] A3727 is propor- 
tional to that of [Olll] A5007 (10% - 30%, Ferland & 
Osterbrock 1986; Ho et al. 1993a, 1993b). Since [Olll] 
A5007 shows a significant BE, [O II] A3727 is unlikely to 
show a radically different trend. So it is safe to conclude 
that the [O II] A3727 and Ha" BE could be partly (but 
not totally) produced by SF contamination. 

A third effect we can see in Fig. 1 is the optical contin- 
uum getting bluer towards higher luminosity. This would 
leverage the continuum and lower the EW if the line flux 
remains unchanged. This could arise if the NLR sees a 
filtered spectral-energy distribution (SED) (Kraemer et 
al. 1998; Groves et al. 2004a, b) so that a change in the 
ionizing spectrum of an AGN would not change the SED 
the NLR sees much. The [Ncv] A3425 , [Oil] A3727 and 
[Nc III] A3870 lines arc most likely to be influenced by a 
continuum shape difference effect because of their shorter 
wavelength. A possible explanation of the bluer color of 
the UV-optical SED is the anisotropy of the continuum 
emission. Because the accretion disc is optically thick, 
we will see a dimmer continuum when viewing it edge-on. 
(This is the combined result of the cos i projection effect 
and the disk equivalent of "limb darkening".) We could 
also preferentially be seeing the inner, high-temperature 
part of the disc when viewing face-on. These effects will 
combine to give a higher observed luminosity with a face- 
on viewing angle. The reddening of spectrum depends on 
viewing angle too (Keel 1980; de Zotti & Gaskell 1985; 
Zhang et al. 2008). Gaskell et al. (2004) used radio 
orientations to get AGN reddening curves and found the 
continuum shape is profoundly affected by reddening for 
all but the bluest AGNs. Because of these effects sources 
with small inclination would both show large luminos- 
ity and have low reddening. An alternative explanation 
of continuum shape difference is host galaxy contami- 
nation. Shen et al. (2011) made composite spectra of 
different ALa(5100) luminosity bins and found that the 
UV part of all the composite spectra are similar while the 
optical part flatten with decreasing luminosity. They in- 
terpreted this trend as due to host galaxy contamination 
in the optical region of the spectrum in low-luminosity 
AGNs. Similar argument is given by Stern & Laor et al. 
(2012a). Even though we have rejected objects with sig- 
nificant stellar light contributions, weak absorption lines 
can be spotted on the final composite spectrum of the 
lowest luminosity bin because of its extremely high sig- 
nal to noise ratio. Our data cannot distinguish different 
mechanisms that give rise to the continuum shape dif- 
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ference, and this is beyond the scope of this paper. Af- 
ter correcting the deviations listed above, our conclusion 
that the slopes of BE for different narrow lines are nearly 
constant is further strengthened. 

4.2. Possible causes of the NLR Baldwin effect 
4.2.1. Softening of the ionizing continuum 

It has been argued on both observational and theo- 
retical grounds that there is a softening of the ionizing 
continuum with increasing luminosity (see, for example, 
Binette et al. 1989; Netzer, Laor, & Gondhalekar 1992; 
Zheng & Malkan 1993; Wang & Lu 1998; Korista et al. 
1998). If this is indeed the case, an important prediction 
of this is the ionization potential - BE slope relationship 
(Korista et al. 1998). It is successful in explaining the 
broad line BE, but a major failing is that it fails to ex- 
plain the constant BE slope of narrow lines. It is already 
known that the NLRs of different AGNs are similar in 
the sense that the line ratios show less than 0.5 dex dif- 
ference from object to object. (Koski 1978; Veilleux & 
Osterbrock 1987; Veilleux 1991a, 1991b, 1991c; Veron- 
Cetty & Veron 2000; Dopita et al. 2002; Gorjian et al. 
2007). Kraemer et al. (2000), Dopita et al. (2002) and 
Groves et al. (2004a, b) proposed a NLR model where 
dust regulates the incident ionizing spectrum so as to 
keep the ionization parameter in the NLR constant. In 
this model, variation of the SED is filtered by dust in the 
NLR and thus the continuum shape seen by the NLR is 
dominated by the effect of the dust rather than by in- 
trinsic changes in the SED. This leads to more stable 
conditions in the NLR. This could be the cause of the 
lack of a dependence of the slope of the BE on the ion- 
ization energy. 

4.2.2. Luminosity-dependence Covering Factor 

A luminosity-dependent covering factor is a natural 
explanation of both a BLR and NLR BE. However, to 
explain the ionization dependence, the variation in cov- 
ering factor must have a different dependence for lines 
of different ionizations. This would not be a surprise for 
the BLR since there is strong radial ionization stratifi- 
cation and the highest ionization lines are an order of 
magnitude closer to the black hole than the lowest ion- 
ization lines (see Gaskell 2009 for a review). Shields, 
Fcrland, & Peterson (1995) suggested that the broad- 
line BE could be caused by a luminosity-dependent cov- 
ering factor for clouds that are optically-thin to pho- 
tons with wavelength less than 912A. The optically-thin 
clouds which have small column densities, preferentially 
emit high-ionization lines. However, comparison of line 
profiles shows that the BLR BE is due to changes in 
the low velocity BLR gas (Francis et al. 1992) rather 
than changes in the very broad component . Furthermore, 
Snedden & Gaskell (2007) argued that optically-thin gas 
does not make a substantial contribution to the BLR. 
Nevertheless, given the strong radial ionization stratifi- 
cation of the BLR, and the decreasing covering factor 
with ionization (e.g., Francis et al. 1992), an additional 
luminosity-dependence of the covering factor could ex- 
plain the BE and why the high-ionization lines should 
show a steeper slope than low-ionization lines. While 
this could be consistent with the BLR, it is not obvious 
how such an explanation can be reconciled with the lack 



of an ionization-dependence we find here for the BE. 

4.2.3. A Disappearing NLR? 

The NLR size is correlated with the luminosity of an 
AGN as Rnlr oc L 5 (Bennert et al. 2002; Schmitt 
et al. 2003b; Bennert et al. 2006a, b; Greene et al. 
2011). So in high-luminosity AGNs, the NLR may pos- 
sibly turn from ionization-bounded to matter-bounded, 
and the luminosity of narrow lines would cease increas- 
ing with AGN luminosity. This model could in principle 
explain part of the BE of NLR (Croom et al. 2002) but 
it would be too rash a conclusion that it is the origin 
of the BE because the bulk of the NLR emission origi- 
nates from within the central few tens or hundreds of pc 
(Schmitt et al. 2003a) so is unlikely to exceed the scale 
of the galactic bulge. Besides, the NLR has no natural 
size because it has no definite edge, so some cut in sur- 
face brightness or line-ratio is needed to define the size 
(Schmitt et al. 2003b; Bennert et al. 2006a, b; Greene et 
al. 2011). This makes the interpretation of size compli- 
cated. Netzer et al. (2004) argue that that Rnlr oc L 5 
is theoretically sound yet this relationship must break 
down for Rnlr exceeding a few kpc. They found that 
high-z AGNs have NLR sizes no larger than lOkpc. It 
is also well known that the the NLR is stratified so that 
the high-ionization lines rise from the inner part of NLR 
while low-ionization lines from further away (Veilleux et 
al. 1995; Bennert et al. 2006a, b). So the "disappear- 
ing NLR effect" , if exists, could preferentially influence 
low-ionization and low-critical density lines. This again 
stands against the constant slope we find. 

4.2.4. Continuum Variation 

Because different lines all seem to share a similar BE 
slope, it is natural to think that it might be the contin- 
uum, rather than the emission lines, which is the cause. 
Jiang et al. (2006) proposed that the continuum varia- 
tion can produce a weak BE (f3 = —0.05 ± 0.05) if the 
light-crossing time for the region emitting narrow FeKa 
exceeds the variability timcscale for the X-ray contin- 
uum and the amplitude of variability anti-correlates with 
the luminosity. Shu et al. (2012) found a strong anti- 
correlation between the EW of the narrow FeKa line and 
L x (EW/(EW) oc (L/(L))-°- S2±0 - w (where () means 
time-averaged values) consistent with the X-ray BE ex- 
pected in an individual AGN if the narrow-line flux re- 
mains constant while the continuum varies. For a NLR 
whose light-crossing time is about 10 3 yr, the narrow-line 
luminosity can safely be assumed to be constant. The 
amplitude of variability of the ionizing continuum can 
be obtained from the monitoring of AGNs after correc- 
tion for the constant host galaxy light contribution. As 
is well known, the amplitude of variability of an AGN in- 
creases with time. This can be most readily seen from the 
structure functions of AGNs (the variance as a function 
of time interval between observations). We are inter- 
ested in variability on the longest timescales (timescales 
similar to the NLR light-crossing time). We do not, of 
course, have monitoring on such long timescales, but we 
can estimate the amplitude from structure functions. Cid 
Fernandez et al. (2000) have presented structure func- 
tions for bright AGNs. These are all consistent with the 
amplitude of variability rising as the time interval in- 
creases to a characteristic time of a year or a few years 
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and then remaining constant. Collier & Peterson (2001) 
got a similar result for lower-luminosity AGNs (but with 
a shorter characteristic time) and also find that the forms 
of the UV structure functions are similar to the optical 
ones. There are indications that the structure function 
increases gradually on much longer timescalcs, but data 
on this part of the structure function are limited. Ob- 
served variability thus gives a lower limit to variability 
on the light-crossing timescale of the NLR. 

In order to minimize the host galaxy light contribu- 
tion one needs to go to as short a wavelength as possible. 
This means using U -band observations or space ultravio- 
let observations. For example, for NGC 4151 Merkulova 
(2006) found an observed peak-to-peak [/-band ampli- 
tude of a factor of 7.5. UV variability can be a lot larger. 
For example, the peak-to-peak amplitude at 1300A for 
Fairall 9 is a factor of 25 (Koratkar & Gaskell 1989). 
Inspection of the results of long-term UV monitoring of 
a number of AGNs with the IUE satellite (Koratkar & 
Gaskell 1989, Koratkar & Gaskell 1991a,b, Clavcl et al. 
1991, and O'Brien et al. 1998) give a median UV peak- 
to-peak variability of a factor of 7, but the upper quar- 
tile is a factor of 13. After correction for host-galaxy 
light, the amplitudes of optical variability are similarly 
large. For example, the NGC 4151 photometry of Lyuty 
& Doroshenko (1999) gives a peak-to-peak amplitude of 
a factor of about 25 to 30. 

Continuum anisotropy (Wang & Lu 1999) can produce 
an effect similar to that caused by the actual variability 
of the continuum. It was recently found by Risaliti et al. 
(2011) that the optically-thick disc emission responsible 
for the continuum and isotropic [O III] emission will pro- 
duce the EW([OlIl]) distribution very well. Variability 
and continuum anisotropy are, in practice, indistinguish- 
able, so we can consider them together. 

We made Monte-Carlo simulations to explore whether 
continuum variability can explain the BE slope. In a 
flux-limited survey, the luminosity will show an approx- 
imately log-normal distribution. We generated an ar- 
tificial sample of sources with a similar distribution of 
ALa(5100) as the observed AGN sample (0.35 dex here). 
We set the EW of the artificial-sources to be the mean 
EW of the whole sample, and we add a Gaussian of 
0.22 dex to the EW to account for the intrinsic disper- 
sion. UV and optical AGN variability is approximately 
log- normal - i.e., it looks normal when plotted in mag- 
nitudes (see Gaskell 2004 for a discussion of log-normal 
variation of AGNs). The continuum variation was simu- 
lated by adding a Gaussian with a = 50% (0.18 dex) to 
the continuum while the emission line flux was kept un- 
changed. This is an over-simplification, but it can help us 
to gain some insight into the effect that continuum vari- 
ability brings. We generated 4000 sources in each round 
and measured the slope of BE using the 4000 sources. 
An example of the simulation is shown in Fig. 4. With a 
pcak-to-peak continuum variation of a factor of three, the 
simulation can produce a BE slope of —0.2 ±0.01 while a 
factor of six variation produces a BE slope of — 0.3±0.01. 
The simulated distribution of EW-ALa(5100) can be seen 
to be very similar to what is observed. Obviously, this 
model will produce a similar slope for every narrow line 
if the continuum change with similar amplitude in the 
wavelength range we concern. It should be noted that 



due to variation in continuum slope, not all lines show 
the same BE slope, but depends on wavelength. We as- 
sume a constant narrow line flux during the continuum 
variability. This is an approximation for the continuum 
variations on time scales of much shorter than the light 
travel time over the NLR because at such short time 
scale, the NLR has little response to the continuum vari- 
ations. For variations on longer time scales, one needs 
to properly convolve the continuum variations with the 
response of emission lines (transfer function). By con- 
sidering the latter response, the variations of equivalent 
width will be somewhat smaller, but this is equivalent to 
requiring a somewhat larger continuum variability am- 
plitude. 

Conversely, since continuum variability or the equiv- 
alent of continuum variability will inevitably produce a 
BE, the BE slope we observe could give an upper-limit 
on the few 10 2 yrs to 10 3 yrs variability of AGN optical 
continuum. If the variability amplitudes exceed the fac- 
tor of three a steeper BE will emerge. This might seem 
at first sight to be at variance with the amplitudes of 
UV variability and structure functions observed, but it 
must be remembered that the equivalent width is the ra- 
tio of the line flux to the observed optical continuum and 
the observed continuum has a substantial starlight con- 
tamination. An example of how the apparent continuum 
variability is substantially less than the real variability is 
shown in Fig. 4 of Gaskell et al (2008). Host galaxy light 
limits the apparent peak-to-peak continuum variability 
to about a factor of three (i.e., a rms variability of a few 
tenths of a magnitude). This is just what is observed for 
the PG AGNs studied by Cid Fernandez et al. Since the 
structure functions are relatively flat after a year or so 
(Cid Fernandez et al. 2000), the variability effect on the 
BE should be apparent in a few years. An obvious test 
of this is to reobserve AGNs after a few years. 

4.3. Drawbacks 

In summary, the combination of a 50% variation of 
ALa(5100) and a log-normal distribution of luminosity 
will naturally produce a -0.20 slope of BE for every nar- 
row line, as we observe. The model we employ is ob- 
viously an oversimplification. To make a more realistic 
model, we need a dedicated treatment of sample selection 
effects and to make a more realistic assumption of the 
variability including amplitude, variation form, and their 
dependence on wavelength etc, about the light curve. De- 
spite these simplifications, our results do show that con- 
tinuum variation will inevitably produce a similar BE for 
each narrow line. Meanwhile, while the model predicts 
the same BE slope for each narrow line, the differences 
in slope are still significant due to the small error bars. 
There must therefore be some additional factor at work. 
It has been shown in Fig. 3 that the differences do not 
correlate with the ionization potential or critical density, 
so there must be other factors taking effect. A deeper 
exploration of other factors like ionization slope, NLR 
geometry and a more realistic model are needed to make 
progress, but these are beyond the scope of this paper. 

5. CONCLUSION 

We have constructed a sample of 3677 z < 0.8 radio- 
quiet AGNs from Sloan Digital Survey Data Release 4 
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that span three orders of magnitude in luminosity to ex- 
plore the relationship between equivalent width of nar- 
row lines and ALa(5100). Wc have computed composite 
spectra for each ALa(5100) bin to enhance the S/N ra- 
tio. Ha b and H(3 b as well as Fell were subtracted to get 
an accurate measurement of narrow emission lines and 
continuum. We find that most narrow lines show a sim- 
ilar BE slope of about -0.2 while the large deviation of 
[Nil] A6583, [Oil] A3727, Ha™ and [NeV] A3425 might 
be explained by a metallicity effect, SF contamination, 
or the UV-optical continuum difference. The slope does 
not show any correlation with ionization energy and crit- 
ical density. We propose that the combination of a 50% 
variation of the continuum near 5100A and a log-normal 
distribution of observed luminosity distribution will nat- 
urally produce a -0.2 slope of BE for every narrow line. 
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TABLE 1 

Ionization energy, Critical density and Slope /3 of log EW - log AL A (5100) Table 



Line 




log ra c 


AL A (5100) slope 


(i) 


(2) 


(3) 


(4) 


[Ncv] A3425 


97.12 


7.3 


-0.31 ± 0.015 


[On] A3727 


13.62 


3.5 


-0.37 ± 0.011 


[Nciii] A3870 


40.96 


5.5 


-0.26 ± 0.020 


H/3" 


13.6 


Infinity 


-0.28 ± 0.015 


[Om] A5007 


35.11 


5.8 


-0.21 ± 0.016 


[O i]A6300 





6.3 


-0.16 ± 0.030 


Ha" 


13.6 


Infinity 


-0.29 ± 0.033 


[Nil] A6583 


14.5 


4.82 


-0.10 ± 0.014 


[Sn] AA6717, 6731 


10.36 


2.3 


-0.20 ± 0.016 



Columns from left to right: (l)The name of the line. (2) The ionization energy needed to create the ion. (3) The 
critical density of specific line. (4) The log EW - log ALa(5100) slope. 
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Fig. 1. — Normalized composite spectra with the broad Ha, the broad H/3 and Fell subtracted are shown for luminosity bins of 
AZogAL A (5100) = 0.3 dex, starting from the top with log AL A (5100)=43.35 [ergs -1 ]. The spectra are normalized to the [4200A,4300A] 
window and shifted vertically to show the weakening of lines with luminosity more clearly. The lines we concern with are marked with 
dashed lines labeled at the top. The embedded plots show the fitting results for the Ha+[N II] + [S II] region. The left embedded panel is a 
composite spectrum for log AL A (5100)G [43.2,43.5] [ergs -1 ] and the right one is for log AL>,(5100)G [44.4,44.7] [ergs -1 ]. 
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Fig. 2. — Line equivalent widths: W\, against 5100A continuum luminosity: AL^(5100). We show weighted linear regressions (green 
lines) for each line and give the BE slopes as well as their errors in the upper right corners in each panel. In the narrow H/3 panel we also 
include the narrow component of Ho (dashed line) for comparison. 
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Fig. 3. — Panel (a) BE slopes of different lines against the ionization energy: Xion needed to create the ion. The diamonds and triangles 
are our data; the red triangles are lines that are contaminated by star formation. The data from Dietrich et al. (2002) and Keremcdjicv 
ct al. (2009) are shown with solid blue rectangles and purple crosses respectively, and the typical error bars are shown in left up corner. 
Panel (b) BE slope against critical density(ra c ) for the same lines in as in panel (a). The recombination lines are marked with right arrows. 
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Fig. 4. — The variability-driven BE simulation result. Panel(a):The simulated BE, the slope of BE is shown on the right up corner. 
Panel(b): The observed EW-log AL A (5100) distribution of [OIII] A5007 from our SDSS DR4 sample. The slope of BE is shown on the 
right up corner too. 



